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Abstract 

We propose here a new model to describe biological invasions in the plane 
when a strong diffusion takes place on a line. We establish the main proper- 
ties of the system, and also derive the asymptotic speed of spreading in the 
direction of the line. For low diffusion, the line has no effect, whereas, past 
a threshold, the line enhances global diffusion in the plane and the propa- 
gation is directed by diffusion on the line. It is shown here that the global 
asymptotic speed of spreading in the plane, in the direction of the line, grows 
as the square root of the diffusion on the line. The model is much relevant 
to account for the effects of fast diffusion lines such as roads on spreading of 
invasive species. 

Keywords: KPP equations, react ion- diffusion system, fast diffusion on a line, 
asymptotic speed of propagation. 
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1 The model 

It has long been known that fast diffusion on roads can have a driving effect on the 
spread of epidemics. A classical example is the spread of the "Black death" plague 
in the middle of the 14th century, considered to be one of the most devastating in 
human history. This pandemics is known to have spread first along the silk road. 
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After reaching the port of Marseilles, carried by merchant boats from Crimea, it 
spread northwards in Europe at a fast pace along the commercial roads connecting 
the cities that had trade fairs. It then also spread more slowly away from the roads, 
inland, bringing about a dramatic invasion. See, for instance, the account by [15]. 

More recently, it has been observed that invasive species such as the Procession- 
ary caterpillar of the pine tree in Europe, have been moving faster than anticipated. 
One plausible explanation is that enough individuals might have been carried on 
further distances than usual by vehicles travelling on roads going through infested 
areas. In the same vein, the invasion of the Aedes albopictus mosquito (also known 
as 'Asian tiger mosquito") is a concern of public health in Europe, see [3]. The 
invasion by this insect is driven by roads. Rivers may accelerate the spread of plant 
pathologies: see for instance [11]. 

Another example of the effect of lines on propagation in open space comes from 
the observation of the population of wolves in the Western Canadian Forest. GPS 
observations reported by McKenzie et al. [14], 2012, suggest that wolves move and 
concentrate on seismic lines. These are straight lines (with a width of about 5m) 
used by oil exploration companies for testing of oil reservoirs. A kinetic model 
leading to an anisotropic diffusion equation was recently proposed by Hillen and 
Painter [9] ; see also [8] . Performing a Chapman- Enskog analysis (with a hyperbolic 
scaling) they derive a drift-diffusion equation where the drift term is directed along 
the lines, and an anisotropic diffusion tensor. They carry out numerical simulations 
that show qualitative agreement with observations: the lines clearly enhance motion. 

From a mathematical point of view, this raises the question of whether the 
inclusion of a line with fast diffusion may affect the overall invasion speed for species 
that usually thrive and reproduce in regions with moderate diffusion capacities. In 
addition, a natural aim is to quantify this effect. The present paper addresses these 
questions and derives precise answers. 

We propose here a model where the two-dimensional environment includes a line 
on which fast diffusion takes place while reproduction and usual diffusion only occur 
outside this line. Thus we consider the line {(x,0) : x G M} in the plane M 2 . For 
the sake of simplicity, we will refer to the plane as "the field" and the line as "the 
road" . But clearly, the model we propose here is a very general one that is relevant 
not only to ecology and population dynamics but also for a wide array of situations 
in biology. For a single species, we consider a system that combines the density of 
this population in the field v(x,y,t) and the density on the line u(x,t). Exchanges 
of populations take place between the road and the field. Namely a fraction v of 
individuals from the field at the road, that is, v(x, 0, t), join the road while a fraction 
[i of the population on the road, u(x,t), goes in the fields. It is assumed that the 
population in the field is subject to a logistic type of growth resulting in a Fisher- 
KPP type of reaction term f(v) in the field. We assume that no such reaction is 
relevant on the road. The diffusion coefficient in the field is represented by d and 
in the road by D. By reasons of symmetry, it will be enough to consider the half 
plane Q := {(x, y) : x G M, y > 0}. Translating into equations all that we have 
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just described, we are led to the following system: 

d t u — Dd xx u = uv(x, 0, t) — [iu x G R, t > 

d t v - dAv = f(v) (x, y) en, t>0 (1.1) 

—dd y v(x, 0, t) = fiu(x, t) — uv(x, 0,t) x G R, i > 0. 

The various parameters are supposed constant D > 0, d, fi, v > and are given. 

System (1.1) incorporates the case where the population also lives in the lower 
half-plane R x (— oo,0) and is symmetric with respect to the x axis. Indeed, since 
u takes the contribution by v from both sides of the road and split its contribution 
to the upper and lower half-plane, the problem reduces to (1.1) with v replaced by 
2v in the first equation and /i by ^ in the last one. This change does not affect the 
linearised system and then, as we will see in the sequel, the spreading properties of 
solutions. Indeed, one come back to the original system by simply multiplying v by 
2. 

Our assumptions on the reaction term / are that / G C 1 ([0, 1]) and satisfies 
/(0) = /(l) = 0, VsG(0,l), < f(s) < f(0)s. 

We extend it to a negative function outside [0, 1]. 

Finally, rescaling the time by a factor 1/r/ in (1.1), we get a multiplicative term 
\jv in front of u, /i, D, d, f. Thus, we assume in the rest of the paper that v — 1; 
this does not entail any loss of generality. 

The quantity 

Ckpp := 2 v / df'(Q) 

is the spreading velocity in the usual KPP equation [1], 

u t - du xx = /(«). 

As we will see, the problem has the structure of a monotone system [10]. See [2], 
[7] for results of front propagation in monotone systems, with which our problem 
shares many features. The originality here lies in the 1D-2D coupling. The main 
result of this paper is the 

Theorem 1.1 (i). Spreading. There is an asymptotic speed of spreading c* = 
c*(n,d,D) > such that the following is true. Let (u,v) be a solution of (1.1) with 
a nonnegative, compactly supported initial datum (uq,vq) ^ (0,0). Then: 

• for all c> c*, we have lim sup (u(x, t), v(x, y, t)) = (0, 0). 

*^ +00 \x\>ct 

• For all c < c* , we have lim inf (u(x, t), v(x, y, t)) = (1///, 1). 

t->+oo |x|<ct 

(ii). The spreading velocity. If d and \i are fixed, and D varies in (0, +oo), the 
following holds true. 

• If D < 2d, then c*(/i,d, D) = ckpp- 

• If D > 2d, then c*{n,,d,D) > ckpp and lim cJfi, d, D) / \/15 exists and is a 
positive real number. 
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Note that, in the statement of Theorem 1.1, the convergence holds pointwise in y. 
A more precise study of the domains where convergence to 1) or (0,0) holds 
is outside the scope of this paper, and will be done elsewhere. 

In the recent years, there have been a considerable number of works on prop- 
agation in heterogeneous media, and how heterogeneities may enhance, or block, 
propagation. See for instance [16], [4]. See also [5] for many other references. 



2 Conservation of total population 

If (u,v) is a solution of (1.1) with / = 0, the quantity ||w(-,t)||ii(R) + ||i>(-, t) ||ii(n) 
does not depend on t. To see this, suppose that u and v decay faster than some 
exponential functions at time t — 0. Anticipating on the next sections, we assert 
that this property still holds for t > 0, also for the derivatives of u and v, owing to 
parabolic estimates. We can therefore integrate by parts the first two equations in 
(1.1) and we find: 

pT p+oo 

||u(-,T)|| L i( R ) - ||n(-,0)|| L i (M) = / / (v(x,0,t) - /iu(x,t))dxdt, 

JO J-oo 

pT p+oo 

\\v(-,T)\\ L i( n) - |K-,0)|| x i ( n) = -d / / d y v(x,0,t)dxdt 

J0 J -oo 
T p+oo 

/ (fiu(x, t) — v (x, 0, t))dx dt. 

J-oo 

Whence, 

\\u(-,T)\\ L i m + |K,T)||£i ( n) = ||n(-,0)|| L i (R) + ||v(-,0)||£i (n) . 



Biological interpretation. Our model is consistent with the conservation of the 
total population in the case of zero natality/mortality rate. The exchange between 
the line and the open plane exactly compensate each other as is natural. The only 
modification in the total size of the population comes from the "effective birth rate" 
represented by the term f(v). □ 



3 The Cauchy problem 

In this section, we derive the existence and uniqueness result for the Cauchy problem 
associated with (1.1). We prove it when Q is the + 1-dimensional half space and, 
for later use, with the addition of a drift term qeR N ,r G R N+1 . Namely, denoting 
the generic point in Q by (x, y) G M> N x (0, +oo), we study the problem 

d t u - DA x u - q ■ V x u = v(x, 0, t) - fiu x G R N , t > 
d t v - dAv -r-Vv = f(v) (x, y) G Q, t>0 (3.1) 

-dd y v(x, 0, t) = im(x, t) - v(x, 0,t) x G R N , t > 0, 
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combined with the initial condition 



u\t-o = uq in M. N 

10 (3.2) 
v | t=0 = v in fi. 

We always assume that Uq and Vq are nonnegative, bounded and continuous. We 
look for solutions satisfying (3.1)-(3.2) in the classical sense. 

The system (3.1) is not very standard because the coupling appears in the Robin 
boundary condition for v. Therefore, well-posedness has to be proved. 

Proposition 3.1 The Cauchy problem (3.1)-(3.2) admits a unique nonnegative, 
bounded solution. 

The proof of the existence part is given in Appendix A. We derive the uniqueness 
part by showing, more in general, that comparison between sub and supersolutions 1 
is preserved during the evolution in Problem (3.1). 

Proposition 3.2 Let (u,v) and (u,v) be respectively a subsolution bounded from 
above and a supersolution bounded from below of (3.1) satisfying u < u and v < v 
at t = 0. Then, either u <u and v < v for all t, or there exists T > such that 
(iLiV) — (u,v) for t < T. 

PROOF. Let / be the Lipschitz constant of /. Define the following functions: 

(u, v) := (u,v)e~ lt , (u,v) := (u,v)e~ lt . 

They are sub and supersolution of the problem 

dfU — DA x u — q ■ W x u + (/i + l)u = v(x, 0, t) xeR N , t>0 
d t v - dAv - r ■ Vv = h(t, v) (x, y) e fi, t > (3.3) 

v(x, 0, t) — dd y v(x, 0, t) = fiu(x, t) x G R N , t > 0, 

where h(t,v) := e~ lt f(ve lt ) — Iv is nonincreasing with respect to v. Let \ : K. — > M. 
be a nonnegative smooth function satisfying 

X (0) = 0, X ' = in [0, 1], lim X (p) = +oo, 

p— >+oo 

((N -l)(d + D) + \q\ + 2\r\)\x'\ + (2d + D)\ X "\ < 1 in E. 
Then, for e > 0, set 

u(x, t) := u(x, t)+e(x(\x\)+t+l), v(x, y, t) := v(x, y, t)+ne(x(\x\)+x(y)+t+l). 
Using the fact that 

\V x (x ^ X (N))I < \x'(\x\)\, \A x (x ^ X (\x\))\ < \ X "(\x\)\ + (N- l)\ X '(\x\)\, 



: A subsolution (resp. supersolution) is a couple satisfying the system (in the classical sense) 
with the = signs replaced by < (resp. >) signs, which is also continuous up to time 0. 
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one readily checks that (u,v) is still a supersolution of (3.3). Moreover, (u,v) is 
strictly above (u, v) at t — 0. Assume by contradiction that this property does not 
hold for all t > 0. Then, 

T := sup{r > : u < u in R N x [0, r], v < v in ft x [0, r]} G [0, +oo). 

It follows that u < u in IR^ x [0, T], v < v in ft x [0, T]. Moreover, the continuity 
of the functions and the fact that u and v tend to +oo as the space variable goes 
to infinity, uniformly in time, implies that T > and either u — u or v — v vanish 
somewhere at time T. If min R ;v(u — u)(-,T) = then, since for x e R N , <t < T, 

d t u — DA x it — q ■ V x u + (// + Z)u > 0, t) > v(x, 0, t), 

the parabolic strong maximum principle yields u = u in M. N x [0, T], which is im- 
possible. Thus, min R jv(n — n)(-,T) > and mm-^(v — v)(-,T) = 0. The strong 
maximum principle implies that the latter cannot be attained inside Q. Thus, nec- 
essarily, (v — 0, T) = 0, for some £ e M. N . This case is ruled out too, because it 
implies 

> (v - 0, T) - d9„(v - 0, T) > y,[u - u)(£, T) > 0. 

We have shown that (u, £>) is above (u, v) for all t > 0, and then n < u and u < U 
due to the arbitrariness of e. Suppose now that there exists a contact point between 
u and u at some time T. Applying the strong maximum principle to the first equa- 
tion in (3.1) we derive (n, v) = (u,v) for t < T. If u < u for all t, then the same 
arguments as before show that v < v. □ 

Biological interpretation. If two pairs of population densities ("road" and 
"field") are initially ordered, the exchanges between them always occur in such 
a way that the order will continue to hold separately for each population at the mi- 
croscopic level. This effect is somewhat more surprising than the mass conservation 
we discussed earlier. Actually, such a monotonicity property is not a priori included 
in the model but rather comes out as an interesting property. □ 

The above comparison principle immediately extends to generalised sub and 
super solutions, given by the supremum of subsolutions and the infimum of super- 
solutions respectively. In Section 6, we will need a comparison principle for a more 
general class of subsolutions, obtained as the supremum of two subsolutions in a 
given set and extended outside to the one which is larger on the boundary. In the 
case of a single parabolic equation, the comparison principle always holds for such 
kind of generalised subsolutions, but this is no longer true for our system (3.1). We 
indeed need an extra assumption. 

Proposition 3.3 Let E C M. N x (0, +oo) and F C Vt x (0, +oo) be two open sets and 
let (ui,Vi), (u 2 ,V2) be two subsolutions of (3.1) bounded from above and satisfying 

U\ < u 2 on {dE) n (R N x (0, +oo)), v x < v 2 on (dF) n (ft x (0, +oo)). 
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If the functions u, v defined by 



u(x, t) : = 



ma,x(ui(x, t), u 2 (x, t)) if (x,t) G E 
U2(x,t) otherwise, 




m&x(v 1 (x,y,t),v 2 (x,y,t)) if(x,y,t) G F 
v 2 (x,y,t) otherwise, 



u(x,t) > u 2 (x,t) =>■ v(x, 0, t) > Vi(x, 0, t), 

v(x, 0, t) > v 2 (x, 0, t) =>- u(x,i) > ui(x,t), 

then, any supersolution (u,v) of (3.1) bounded from below and such that u < u and 
v<v at t = 0, satisfies u<u and v<v for all t > 0. 

Proof. Since u > u> u 2 and v > v > v 2 at t — 0, Proposition 3.2 yields u > u 2 
and v > v 2 for all t. Assume by contradiction that (u, v) is below (u, v) somewhere. 
Using the same transformations as in the proof of Proposition 3.2, with e small 
enough, one reduces to the case where {u\, v\), {u 2 , v 2 ) are subsolutions of (3.3) and 
(u, v) is a supersolution of (3.3) which is strictly above (u 2 ,v 2 ) for all t, is strictly 
above (u, v_) for t less than some T > and touches it somewhere at t — T. We 
argue differently depending on the fact that u touches u or not at t — T. 

Case 1. u(i, T) = T) for some ^ e R N . 
Since u 2 {^,T) < u(^,T) = m(^,T) = ui(£,,T), it follows that (£,T) G E and there 
exists 5 > such that u 2 < U\ = u in Bg(£) X (T — 5, T + 5) C -E. Let us call 
Q := -£>£(£) x (T — 5, T). For (x, t) G Q, the hypothesis on (m, v) yields Vi(x,0,t) < 
v(x,0,t) <v(x,0,i). Whence, 

d t ui — DA x Ui — q ■ V x ui + (/i + l)u\ < vi(x, 0, t) < v(x, 0, t) in Q. 

But «i < m on the parabolic boundary of Q and then the strong maximum principle 
yields U\ =um.Q, which is impossible. 

Case 2. u<u for t = T and v(£, rj, T) = v(^, 77, T) for some (£, 77) G £1. 



Since v 2 (£,r],T) < v(£,rj,T) = v(£,r},T) = ?Ji(^,?7,T), it follows that {£,r},T) G F 



and there exists 5 > such that t> 2 < ^1 in (F>s(!;,ri) HQ) x (T — 5, T + 5). By 
the hypothesis on V\ and t> 2 , this set cannot intersect OF and then it is contained 
in F, because (£,r],T) G F. Whence, we find that v 2 < v x = v < v in Q := 
(Bs(£,, t]) nil) x (T — 5, T). If 7] > 0, the parabolic strong maximum principle implies 
that V\ = v in Q, which is impossible. Therefore, rj = and then the hypothesis on 
v yields «i(^,T) < u((,,T). Consequently, 



0>(v- <)(£, 0, T) - dd y (v - v 1 )(C, 0, T) = /i(n - T) > fi(u - u)(£, T). 



We have reached a contradiction, because we are in the case u<u for t = T. □ 

Remark 3.4 It is clear that in the statement of Lemma 6.2 one can require that 
(u\, v 1) is a subsolution of the three equations of (3.1) only in E, F and for {x, 0, t) G 
dF respectively. 
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4 Long time behaviour 



We derive a Liouville-type result for stationary solutions of system (3.1) in the 
absence of drift, that is 

-DA X U = V(x, 0) - fill xeR N 
-dAV = f(V) (x,y)eQ (4.1) 

-dd y V(x, 0) = fj,U{x) - V(x, 0) xe R N . 

Proposition 4.1 The unique nonnegative, bounded solutions of (4.1) are (U, V) = 
(0,0) and (U,V) = (l//i,l). 

Proof. Let (U, V) be a nonnegative, bounded solution of (4.1). By the last equa- 
tion of (4.1), the statement is proved if we show that V = or 1. Hence, by the 
strong maximum principle, it is sufficient to show that if V > in Q then V = 1. 
We do it in three steps. 
Step 1. The function V satisfies 

Vr > 0, inf V > 0. (4.2) 

This step is the only part of the proof in which the absence of drift in the system 
is required. Take R > large enough in such a way that the principal eigenvalue of 
—A in Br C K +1 , under Dirichlet boundary condition, is less than f'(0)/d. Let 
tp be the associated principal eigenfunction. It follows that there exists 6$ > such 
that for e < Eq, —dA(e<p) < f(sip) in Br. As a consequence, for given x G M. , 
y > R, the elliptic maximum principle yields 

V(x, y) <G B R (x, y), V(x, y) > e (p(x -x,y-y). 

This proves (4.2). 
Step 2. V > 1. 

Set m := miV and consider ((x n ,y n )) n such that V(x n ,y n ) —> m. We assume by 
way of contradiction that m < 1 and we distinguish two cases. 

Case 1: (y n )nm tends to 0. 
We set 

U n (x) :=U(x + x n ), V n (x,y) := V(x + x n ,y). 

By standard elliptic estimates, we have that {U n ) n and (V n ) n converge (up to subse- 
quences) locally uniformly to some functions U, V satisfying the same system (4.1). 
Furthermore, y(0,0) = m = min^-V^. Since AV^ = —f(V)/d < in the intersec- 
tion of Q with a neighbourhood of (0, 0), the strong maximum principle and Hopf's 
lemma imply that either V = m oi V > m m. VL and —d y V(0, 0) < 0. The first 
situation cannot occur because f{m) > by (4.2). In the second situation, using 
(4.1) we derive fiU{0) — m < 0. In particular, 

77? 

inf U < inf U < —. 
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Considering now the limit U of translations of U by (a subsequence of) one of its 
minimizing sequences, we reduce to the case where 

U(0) = min U = inf U, —DA X U >m-fiU. 

Hence, evaluating the latter at x — 0, we get the following contradiction: 

< DA x U{0) < //inf U - m < 0. 

Case 2: there exists a subsequence (y nk )ken of (y n )neN with distance r > from 

0. 

It follows from (4.2) that m > 0. The sequence of the translated functions Vk(x, y) : = 
V{x + x nk ,y + y Uk ) converges (up to subsequences) locally uniformly to a function 
V satisfying 

< in B p , y(0,0) = m = minV, 

for some < p < r. This is impossible. 
Step 3. V <1. 

Assume by contradiction that M := sup V > 1. Analogously to the previous step, we 
can reduce to the case where V(0, 0) = maxV = M > 1 and either (U, V) satisfies 
(4.1) or -AV = f(V) in a neighbourhood of (0,0). Since f(M) < 0, the latter 
case is impossible and, by Hopf's lemma, we necessarily have that —d y V(0, 0) > 0. 
Whence (4.1) yields 

M 

supf/ > U(0) > — . 

By usual arguments, the limit U of translations of U by (a subsequence of) one of 
its maximizing sequence attains its maximum sup U at and satisfies 

DA x U(0) > [i sup U — M > 0. 

This is a contradiction. □ 



Theorem 4.1 Let (u,v) be a solution of (1.1) with a nonnegative, not identically 
equal to zero, bounded initial datum. Then, 

lim (u(x,t),v(x,y,t)) ->■ (l//x,l), 

t— >+oo 

locally uniformly in (x,y) G f2. 

PROOF. Take i? > large enough in such a way that the Dirichlet principal 
eigenvalue of —A in Br C 1R 2 is less than f'(0)/d. Call </? the associated prin- 
cipal eigenf unction. Hence, for e > small enough, the function V_ defined by 
V(x,y) := e<p(x,y- R - 1) satisfies V < 1 and -dAV; < f(V) in S R (0,i? + 1). 
Extending V_ by outside B R (0,R + 1), we have that (0, VQ is a generalised sub- 
solution of (1.1) in the sense of Proposition 3.3. It is then easy to check that the 
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solution (u,v) of (1.1) with initial datum (0,V_) is nondecreasing in time. On the 
other hand, the pair (U, V) defined by 

U := max ( sup u , — sup v , — J , V := max ( \i sup u , sup v Q , 1 

\ t si \k n 

is a supersolution of (1.1). It is now sufficient to use Proposition 3.2 to infer that 
the solution (u,v) with initial datum (U,V) is nonincreasing in time. Therefore, 
owing to Proposition 4.1, as t — > +oo, they both converge to the unique positive, 
bounded solution of (4.1): (U, V) = (U). This convergence is locally uniform in 
(x,y) by parabolic estimates. Comparing (u,v) with (u,v), we derive 

]imswp(u(x,t),v(x,y,t)) < (l//i, 1), 

t— >+oo 

locally uniformly in (x, y) G Q. Furthermore, by Proposition 3.2 we know that v > 
in Q x R + . Hence, up to choosing a smaller e in the definition of V_ if need be, it is 
not restrictive to assume that V_(x, y) < v(x,y, 1) for (x,y) G Q. Thus, 

liminf (u(x, t + 1), v(x, y,t + 1)) > lim (u(x, t),v(x, y, £)) = (l//i, 1), 
locally uniformly in (x, y) G Vt. The proof is thereby achieved. □ 



Biological interpretation. Even though there is no birth rate on the "road", the 
presence of a logistic term in the "field" is sufficient to ensure global invasion (road 
and field) by the species modelled by (1.1). This is of course due to the exchange 
term between the field and the road. □ 



5 Exponential solutions 

The upper bound for the spreading speed will follow from the existence of exponen- 
tial solutions for the linearised problem 

d t u — Dd xx u = v(x, 0, t) — iiu x G R, t G E 

d t v - dAv = f'(0)v (x,y)ett,teR (5.1) 

—dd y v(x, 0, t) = pLu(x, t) — v(x, 0, t) x G R, t G R, 

because, by the KPP hypothesis, solutions to (5.1) provide entire super-solutions 
to (1.1). The lower bound will follow from the existence of solutions with complex 
exponent. So, we will be looking for solutions of the form 

(u(t, x), v{t, x, y)) = ( e «(*+rf) ? T e«(*+<*)-/»!/) (5.2) 

where a and 7 are positive constants and (3 is a real (not necessarily positive) 
constant. In this section, we analyse the condition on (a, (3, 7) so that u and v given 
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by (5.2) solve (5.1). The system on (a, /?, 7) reads 

-Da 2 + ca = 7 — fi 

-da 2 + ca = f'(0) + d(3 2 (5.3) 
d/?7 = n — 7 

The third equation of (5.3) gives 7 = /i/(l + Hence, the first equation, with 
unknown a and parameters D, d, (3, /i, has the roots 



2D 



dp )' 



For a^(c, 0) to be real, the parameter (3 has to be chosen in (— 00,— 1/d) U 
[(3d(c), +00) with 



/3d (c) 



c 2 



d(c 2 + 4/iD) ' 

Note that if /3 < — 1/d we have 7 < 0. This case has therefore to be discarded. Also 
note that lim (3d(c) = —1/d. Denote T^ D := {a = a^(c, (3), (3 > (3n(c)}; in the 

(f3, a) plane, T c D := r+ D ur~ D is a smooth curve with leftmost point (/3d (c), c/2D). 
See Fig. 1. 

Let us now deal with the second equation in (5.3); in the ((3, a) plane, it is the 
equation of the circle T Cjd of centre (0, c/2d) and radius 



/3kpp(c) 



2 

KPP 



2d 

In the whole discussion below, it will be necessary to require 

c > ckpp = 2^dfW)- 

This will be assumed without any further mention. The circle T c ^ is the union of 
the two half circles T^ d = {a = a^(c, /?)}, where 

, _ c±^c 2 -c 2 KPP -4d 2 (3 2 
d{ ,P) 2d 

The quantity 

"d(c,0) = — - !3 K pp(c) 

is the exponent at which a travelling wave of (1) will decay in the x variable. So, 
(5.3) amounts to finding c such that T cP) and T c ^ intersect. The following properties 
of a^ and a^ are easily checked: 

• The function satisfies d^a^ < and dppa^ > 0; moreover a^(/3,c) < if 
(3 > 0, and < if (3 e [/3 D (c),0). 

• The function a J satisfies d c a~^ > 0, dpa^ > wherever it is defined; moreover, 
d P pa%(c, (3) < 0. 
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• The function a d is even; moreover we have d c a d < 0, df3pa d > 0, and dpa d > 
for (3 > 0. 

• The function a d is even; moreover we have d c a d > 0, dppa^ < 0, and dpa d < 
for > 0. 

Case 1. D > 2d. For c very large we have 

So, there is c* > ckpp such that T Ctt p, and r Cjt)C ; are tangent. The tangency point 
belongs to T~ d ; moreover, the properties d c a^ > and d c a d < imply that Y c D 
and r c (i intersect as soon as c > c*. Due to the convexity properties of T CjD and r c d , 
these two curves intersect at exactly two points denoted by (/3±, a±) with (3_ < f3 + , 
a_ < a + . The situation is summarized in Figure 1. We can be a little more explicit 
and tell to what part of T C} d the points belong. The point here is whether T^ D and 
T^ d intersect, and the condition is 

D<2d + 5 with 5 := ( max— - -. ) And 2 . 

Indeed = a d admits a solution iff aj(c, (3kpp(c)) > a d (c, (3kpp(c)). The latter 

And/3 KPP (c) > D-2d 



c 2 (l + d/3 KPP (c)) ~ d? 



whence, calling t := 2d(3KPp(c) = \J c 2 — c 2 KPP) we get 

\\id 2 t 



D-2d< 



(t 2 + c 2 KPP )(t + 2Y 

which admits solutions t > iff D < 2d + 5. Actually, if D < 2d + 5, since 

(D - 2d)(t 2 + c 2 KPP )(t + 2) = AfxdH 

has two positive solutions (coinciding if D — 2d + 8), it follows that there exist 
c* < c\ < c 2 such that, for all Ci < c < c 2 , the equation aj(c, (3) = a d (c, (3) has 
a unique positive solution /3. We have that 5 < fid/f'(0) , that D h-> ci(D) is 
increasing, D i-> C2(-D) is decreasing and 

lim cAD) = 0, lim c 2 (-D) = +oo. 

D^2d- D^2d- 



Case 2. < D = 2d. The centre of Y Cy d sits on T Ct c, exactly at the point (0,c/D). 
So, for all c > Ckpp, r cD and T cd intersect exactly twice. In this case, we set 
c* := ckpp- 

Case 3. < D < 2d. Similarly to the discussion of Case 1, there is a largest 
d > ckpp such that, for all c < c', T Ct d and Y c ^ do not intersect. The situation 
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Figure 1: Case D > 2d; c < c* (left), c = c* (middle), c > c* (right). 



here is more intricate than in Case 1, and also less interesting as far as spreading 
is concerned: we will indeed see that here, spreading occurs at velocity ckpp- This 
study is therefore left for future investigation. Setting c* := ckpp, we see that for 
c > ckpp, the centre of T c ^ lies on the right of r C) £>. This shows that the exponential 
function associated with a = c/2d, = is a supersolution of (5.1). 
Cases 2 and 3 are illustrated in Figure 2. 



Biological interpretation. There is a threshold value for D here, at D — 2d. 
At first sight, this might appear somewhat surprising as one might have expected 
for instance D = d to be the threshold (identical diffusion on the road and on the 
field). It turns out that this threshold value is directly related to the fact that in 
our model there is no reaction on the road. As a matter of fact we may consider a 
related model where one allows reproduction on the road with the same rate as in 
the field. In this case, the new system reads: 




A 



A 



Figure 2: Cases D = 2d (left) and D < 2d (right). 




(5.4) 
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Adapting the arguments of the next section one can show that then, the threshold 
occurs at D = d. A detailed study of the system with a reaction term on the road 
too will be done elsewhere. □ 



6 Asymptotic spreading 



Because of the comparison principle (Proposition 3.2), the standard way to prove 
that solutions spread at least at speed c* is to find a compactly supported generalised 
stationary subsolution in the moving framework at velocity c < c*. Here, generalised 
subsolution is in the sense of Proposition 3.3. We consider the linearised system 
penalised by 5 > in the moving framework: 



Dd xx u + cd x u = v(x, 0,t) — jiu x6l, teM 

{x,y) en, te 



d t u - 

d t v - dAv + cd x v = (f'(0) - S)v 
-dd y v(x, 0, t) = fiu(x) — v(x, 0, t) 



(6.1) 



x G 



t e 



For D > 2d, the main lemma is the following. 



Lemma 6.1 Assume D > 2d and let c* be as in Section 5. For c < c* close enough 
to c*, there exists 5 > such that (6.1) admits a nonnegative, compactly supported, 
generalised stationary subsolution (u,v_) ^ (0,0). 

Proof. In order to keep the notation as light as possible, we will carry out the 
discussion with f'(0) instead of /'(0) — 5 in (6.1), because all the results will perturb 
for small 5 > 0. 

The first step is to devise a stationary solution of (6.1), not in Q anymore, but 
in the horizontal strip Q L :=Rx (0, L) with L > 0. So, we are solving 



-DU" + cU' = V(x, 0) - fiU 
-dAV + cd x V = f'(0)V 
-dd y V(x, 0) = fiU(x) - V(x, 0) 
{V(x,L) = 



x e R 
(x,y) e tt L 
x e R 
x e R. 



(6.2) 



In a similar fashion as in the preceding section, we look for exponential solutions of 
the form (1, 7(y))e a:c . An easy computation shows that such a solution exists if and 
only if the system with unknowns a and (3 (look for 7 under the form 7ie _/3y + 72e /3y ) 



-Da 2 + ca + 



1 + e-W^dPfji 



e -2 P L + ( X + e -2/3L) d/ g 
-d(a 2 + (3 2 ) + ca 
dp fri - 72) 

-PL + 72e /3L 



7ie 





f(0) 

V- (7i 




72) 



(6.3) 



has a solution. The first equation of (6.3) defines, in the upper part of the plane 
((3, a), a curve T L D , which is symmetric with respect to the vertical axis. Moreover, 
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the set T^ D D {a > 0, > 0} tends to T C)D n {a > 0, > 0} as L ->■ oo. In fact, it 
is the graph a = a^' L (c, /3) with 



<*d L (c,P) 




4(1 + e-^ L )fidD(3 
■ e-^ L + (1 + e- 2 P L )d(3 



The curve is strictly above r cD ; moreover, as L — > oo, converges to a J, 
together with all its derivatives, in every compact subset in (c, /3) avoiding {0 < 0}. 
Thus, due to the implicit function theorem, for large L, the picture is analogous to 
the case 1 of the previous section: there exists a unique G (c^pp, c») such that 
^cD H {/3 > 0} intersects r c(i twice if c > once if c = and never if c < . 
Moreover we have lim = c*. 

To obtain the compactly supported subsolution, fix L > suitably large and set 

h L (c,0) :=a+' L (c,0)-a^(c,0). 

Call (0^ , a*) the tangent point between D D {/3 > 0} and r c z )rf . We have 

^(c^, = 0, -2a := 0^/i L (cf, & L ) < 0. 

Also, set e := d c h L (c^,0^) > 0, 6 := d c ^h L (c^, 0^), and, because we are working in 
a vicinity of (t£,#f): 

e := - c, r := /3 - 
The equation /i L (c, /3) = becomes, for (c, /3) in a neighbourhood of (c%, 0*): 

ar 2 + ^r + ee = <^(r,0 (6.4) 

where tp is analytic in r in a neighbourhood of 0, vanishing at (0, 0) like |t| 3 + £ 2 . For 
small £ > 0, the trinomial ar 2 + d£r + e£ has two roots t± = dbi-y/' (e/a)£ + 0(£); by 
(an adaptation of) Rouche's theorem (see Appendix B), equation (6.4) has two roots, 
conjugate up to the order £, still called r±. They also satisfy r± = ±2a/ (e/a)£+0(£). 
Reverting to the full notations, we see that, for c strictly less than and suitably close 
to , equation (6.3) has a solution (/3, a) with the following properties: has real 
part £-close to 0^, hence positive, and it also has nonzero imaginary part of order 
£ x / 2 . It follows from the second equation of (6.3) that a has nonzero imaginary part 
too. Therefore, one can write 

(0, a) = ((3i + i0 2 , ati + ia 2 ), 

with f3i > and a 2 , 02 ^ 0. We thus obtain a solution (U(x),V(x,y)) = 
(l,^(y))e^ ai+ia2 > to (6.2), with 

nf(v) = ^ (e - ^ - e^ -2 ^) = - 

W h 7i 1 _ e -2 f }L + d/3 ( 1 + 2e -2pLy 

We have arg(7i) = 0(y/£). Furthermore: 
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• Re(£7) > on (— 7r/2a 2 , 7r/2a 2 ) and vanishes at the ends; 

• Re(V) > if and only if 

cos(a 2 x + arg(7i) - f3 2 y) > e _2/3l(i_2/) cos(a 2 x + arg(7i) - (3 2 (2L - y)). 
That is, 

sin(2/3 2 (-^ — 2/)) 

cos(a 2 x + arg(7i) - p 2 y) > &2 ^ L _ y) _ CQS ( 2/ g 2 (L _ g sin(a 2 x + arg( 7l ) - fay). 

Thus, the set where Re(V) > is periodic in the direction e±, with period 
Its connected components intersecting the strip Rx (0, L) are bounded. Since (3 2 = 
by decreasing £ we can make one of them, denoted by F, satisfy the property 
that {a 2 x + arg(7i) : (x, 0) G F} is arbitrarily close to the set [— 7r/2, 7r/2]. Since 
arg(7i) = O(v^), we have that {a 2 x : (x,0) G F} is close to [— 7r/2,7r/2] for £ 
small. We define the following functions: 



v(a:,2/) := 



max(Re(?7(x)),0) if |x| < ^ 
otherwise, 



max(Re(y(x, y)), 0) if (x, y) e F 
otherwise. 



The choice of F implies that the couple (u,v) is a generalised subsolution to (6.1) 
in the sense of Proposition 3.3. □ 

Turn to the case < D < 2d. The main lemma is here the following: 

Lemma 6.2 Assume < D < 2d. For —ckpp < c < ckpp, the conclusion of 
Lemma 6.1 holds. 

Proof. Let |c| < c K pp- For < 5 < /'(0) - c 2 /4d, the equation 

-dAV + cd x V = (f(0) - 5/2)V, (x, y) g M 2 (6.5) 

has a compactly supported subsolution. The construction is classical: setting u := 

we see that the Dirichlet problem 

-d(f>" + ctf>' = (f(0) - 5)<p in (-7r/2w, 7r/2w), 0(±tt/2u;) = 

admits the positive solution (j>($) = e^ c ' 2d ^ x cos(ux). Then, let ipp(y) be the first 
eigenfunction of — d yy in (-R, R); choosing R such that the associated eigenvalue is 
S/2d, it is readily seen that V(x, y) := <f>(x)if}R(y — R— 1) is a solution to (6.5) in the 
rectangle (~tt/2uj, tt/2u) x (1, 2R + 1), vanishing on the boundary. Extending it by 
outside, we find that (0, V) is a generalized subsolution to (6.1) (cf. Remark 3.4). □ 

Proof of Theorem 1.1: the 'spreading' part. Let c* be as in Section 5 and 
let ( e a ( I + c « t ) ) -ye Q ' I+c *^~^) be the exponential entire supersolution of (5.1) - and 
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then of (1.1) - constructed there. By symmetry, (e a<y x+c * t ) ^e a( - x+c **) 0f) is still a 
supersolution. The comparison principle - Proposition 3.2 - then yields 

Vc > c*, lim sup (it(x, t), u (x, y, t)) = (0, 0). 

*^+°°|a:|>ci 

Fix now < c < c* (close enough to c* if -D > 2d) and consider the pair (u,v) 
given by Lemmas 6.1, 6.2. There exists 7o > such that, for < 7 < 70, j(u,v) is 
a subsolution of the problem in the moving framework: 

d t u — Dd xx u + cd x u = v(x, 0, t) — fiu x £ R, £ > 
fifcv - c/At; + cS^u = /(u) (x, y) EVl, t>0 (6.6) 

—dd y v(x, 0, i) = /iu(x, t) — v(x, 0,t) x £ E, £ > 0. 

For < 7 < 70, let (it 7 , u 7 ) be the solution of (6.6) with initial datum j(u, v). Using 
the comparison principle for generalised subsolutions - Proposition 3.3 - we see that 
(tt 7 , u 7 ) is nondecreasing in £. Moreover, since j(u, v) is not a solution of (6.6), it 
follows from the strong comparison principle given by Proposition 3.2 that (w 7 ,t> 7 ) 
is strictly increasing in £. Thus, as t +00, (u y ,Vj) converges locally uniformly 
to a stationary solution (U y , V^) of (6.6) which is strictly larger than j(u,v). There 
exists then k > such that (U y , V 7 ) is above the translated by any h £ (— k, k) in 
the x-direction of 7(7/, v_). By comparison with the translated by h of (w 7 ,f 7 ), we 
infer that (U 7 , V y ) is above the translated by h of itself, that is, it does not depend 
on x. Therefore, U 1 is constant and V y = V y (y) satisfies 

Yy(0)=nU 7 , -dV^ = f(V,) fory>0, tfj(0) = 0. 

This easily implies V 7 = 1 and U 7 = We will now conclude by comparing the 
solution (u, v) with a combination of the functions 

(u 7 (x + ct, t),v 1 (x + ct, y, £)), (m 7 (— x + ct, t),v 7 (—x + ct, y, £)), 

which are also solutions of (1.1). Notice that, up to waiting until time 1, it is not 
restrictive to assume that, at initial time, (u, v) is positive and then it is above 
■y(u, v) for some 7 < 70 . Take T > 2 and < £ < c(T — 2). Let r £ [1, T/2] be such 
that £ = c(T — 2r). By comparison, we get 

V(x, y) £ fi, < £ < r, (w(x, £), v(x, y, £)) > (u 7 (x + ct, £), v 7 (x + ct, y, t)), 

whence, in particular, 

(u(x, r), v(x, y, r)) > (u 7 (x + cr, r), w 7 (x + ct, jy, r)) > (u 7 (x + cr, 1), v^x + cr, y, 1)). 

Taking < 7' < 70 small enough (independently on T) in such a way that 
■y'(u(—x),v(—x,y)) < (u 7 (x, 1), w 7 (x, y, 1)), we find that, at t — r, (u,v) is larger 
than the pair 

(ity (— x + c(£ — 2r), £ — r), ty (— x + c(£ — 2r),y, t — r)). 
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Applying once again the comparison principle we get, for given y > 0, 
(u{£,T),v(£,y,T))> (u Y (0,T-r),v y (0,y,T-r)) > (u y (0,T/2),v y (0,y,T/2)). 
From this, we eventually infer 

The reverse inequality follows from Theorem 4.1. This concludes the proof, as the 
negative values of £ can be handled by reflection with respect to the y axis. □ 

Biological interpretation. The result we have derived in this section is a remark- 
able effect of this model. Essentially, it says that propagation on the road can be 
accelerated by a large diffusion, provided that there is an exchange between the 
road and an outside region where reproduction and some diffusion occurs. For lower 
values of the diffusion on the road, namely, D < 2d, the propagation on the road is 
purely driven by that in the field. The invasion velocity then is the standard KPP 
invasion speed. However, when the critical value D = 2d is crossed, a transition to 
a new regime takes place. Namely, the combination of reaction in the field and fast 
diffusion on the line yields an invasion speed strictly larger than the KPP invasion 
speed. One may view this regime as driven by the diffusion on the line. Note that 
it is rather remarkable that even though no reproduction takes place there (lead- 
ing to dispersion if it were an isolated line), the large diffusion drives the overall 
phenomenon. 

The precise value of the threshold diffusion here is D — 2d. In the previous 
section, we saw that when the reaction (effective reproduction) also takes place on 
the road, one is led to system (5.4). For that system, the threshold value is D = d. 
It is interesting to note that while in both cases a large diffusion on the road always 
enhances the spreading speed, when there is pure diffusion and no reproduction on 
the road, this effect starts to be felt at higher values of the diffusion when compared 
to the model in which the reaction holds everywhere. □ 



7 The large diffusion limit of c* 

In this section we investigate the behaviour of the critical speed c*, introduced in 
Section 5, as D goes to oo. Hence, c* is now understood as a function of the variable 
D. We will first show that 

2 2 

v/4/i 2 + (J' (0)) 2 - 2/i < liminf ^ < limsup^ < /' (0). (7.1) 

D^oo D D ^oo D 

The proof of the last statement of Theorem 1.1 will follow easily. 

For D > 2d, we know that c* > ckpp = 2^df'(0). Moreover, it follows from 
geometrical considerations (see figure 1) that 

% < Yd ~~ P Rpp ( c * ) < + v /c * + 4/ijD)) 
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where, we recall, (3kpp(c) = — — 2d KPP is the radius of the circle T Cjd . Whence, 



1 




From the second inequality above we infer that, as D — > oo, c* tends to oo and then 
1 1 




Calling r := the first inequality above yields r 2 < /'(0) + o(l), and then, by 
the second one, 



2/(0) + o(l) < r 2 + r^/r 2 + 4/x < /'(0) + r^r 2 + 4/x + o(l). 
Property (7.1) then follows. 

Proof of Theorem 1.1: the limit of c*. Revert to the algebraic system (5.3). 
Estimate (7.1) enables us to rescale the unknowns c and a as 

c = VDc, a = a/vD, 

with c bounded from below away from and from above independently on D. It 
then follows from the first equation of (5.3) that a is bounded from above too. Drop 
the tildes. Because c and a are now to be found bounded, we drop the term —da 2 /D 
produced by the normalisation in (5.3), and the limiting system to solve is 

, d(3u 
-cr + ca = — — 

1 1 + dfi (7.2) 

a = -(/'(0) + d/3 2 ). 
c 

The last equation of (7.2) is that of a parabola in the (/3, a) plane, whose bottom is 
the point (0, /'(0)/c). Call it - just to be consistent with the notations of Section 
5. The first equation just represents the set Y c D for D = 1. When c is close to 0, 
the parabola r^° d is above this set - which lies below the line {a = 2-y/Jl}. When c is 
close to infinity, because of their respective convexity properties, the curves intersect 
at two points (see Figure 3 below). Thus (same argument as in Section 5) there is 
c* > such that (7.2) has exactly one solution. This is the sought for limit. □ 

Biological interpretation. When the diffusion on the road is very large, the 
spreading along the road also becomes very large. The invasion speed can be made 
as large as one wishes provided the diffusion on the single line is taken sufficiently 
large. This model therefore accounts for the enhancement of invasion speeds in the 
presence of lines with large diffusion coefficients even though no growth of popu- 
lation takes place there. The precise quantification we derive here states that the 
invasion speed grows like the square root of the pure diffusion on the road. It seems 
interesting to go further into this analysis and determine how the prefactor of \f~D 
depends on the various parameters. □ 
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Figure 3: c close to (left), c = c* (middle), c close to oo (right). 

8 Conclusion 

We have introduced a model that describes the effect of a line on which there is fast 
diffusion on the overall propagation of a species that diffuses with another constant 
and reproduces outside this line in a two dimensional framework. 

We have found that this model conserves the population in absence of reproduc- 
tion and mortality and preserves order. Then, we have shown that owing to the 
exchanges taking place between the line and the plane, there is an asymptotic speed 
of spreading which is the invasion velocity along the line. 

We have computed the global asymptotic speed of spreading along the line. 
This is achieved with exponential solutions of the linearised system and compactly 
supported sub-solutions. The asymptotic speed is derived from an algebraic system. 

When D, the diffusion on the road, is less than or equal to 2d, where d is the 
diffusion in the field, there is no effect at all of the road: the propagation takes place 
at the classical KPP invasion speed. In contradistinction with this case, when D is 
larger than 2d, there is an enhancement effect of the diffusion on the road leading 
to a speed higher than KPP. Lastly, this invasion speed is shown to behave like y/~D 
for large values of D. 

Appendix 

A Existence result for the Cauchy problem 

Proof of the existence part of Proposition 3.1. We prove the result for an 
initial datum (uq, v ) which is locally Holder continuous, together with its derivatives 
up to order 2, and satisfies the compatibility condition 

-dd y v (x, y, 0) = iiu (x, y) - v (x, y, 0). 

The regularity of the initial datum is therefore inherited by the solution of the 
Cauchy problem for all time t > 0. The case of a merely continuous initial datum 
can then be handled by a standard regularization technique (see, e.g., [13]). 
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We will obtain a solution to (3.1)-(3.2) as the limit of a subsequence of solutions 
{{u m v n )) n of the following problems: 



pN 

TT"W A 

u n \ t =o = u m 



d t u n - DA x u n - q ■ V x u n + fiu n = v n -i(x, 0, t) x e R N , t > 

N IA.1J 



<9 t i> n - dAt?„ - r • Vw n = / (v n ) (x, y) e fi, t > 

w n (x, 0, t) - dd y v n (x, 0, t) = pu n (x, t) x e R N , t > (A.2) 
v n \ t =o = v in fi, 

starting from Vq. 

Step 1. Solvability of (A.l), (A.2). 

We say that a function i) has admissible growth in z if it satisfies t)\ < 
/3e ff l^ 2 ,for some a, /3 > 0. It is well known that the linear Cauchy problem is uniquely 
solvable in the class of functions with admissible growth in the space variable. If 
v n -i is a continuous function with admissible growth, then problem (A.l) admits a 
unique classical solution u n with admissible growth. In order to solve (A.2), notice 
that it can be reduced to a homogeneous system by replacing v n with v n — vq — 
/j,(u n — Mo). It then follows from the standard parabolic theory that it admits a 
unique classical solution with admissible growth. Let ((u n ,v n )) n denote the family 
of solutions constructed in this way, starting from Vq. 
Step 2. L°° estimates. 
We show, with a recursive argument, that 

Wn G N, 0<u n <—H, 0<v n <H, with H :— max (ii||«o||oo, lkolloo, 1) ■ 

The property trivially holds for n — 0. Assume that it holds for some value n — 1. 
Since and —H are respectively a sub and a supersolution of (A.l), the comparison 

principle yields < u n < -H. It then follows that if is a supersolution of (A.2), 
whence < v n < H. 
Step 3. Wp' 1 estimates. 

By step 1 we know that < t> n _i < H. Thus, applying the local boundary estimates 
to (A.l) we infer that, for any given p, T > and 1 < p < oo, 

IWIvVp'^Bp+ixCO.T)) — CH, 

where B p denotes the N- dimensional ball of radius p and centre and C is a constant 
only depending on N, D, q, p, p, T, p and ||wo||w^(s p +2) (and not on n). Set 
Q p := Bp x (0,p). Since < v n < H too, the estimates yield 



I) < C" (i + IK|U°°(Q p+1 x(0,T)) + ||w n || w 

< C"(l + # + C#), 



with C" only depending on N, d, r, p, p, T, p, U/H^, ||u || w£(b p+1 ) and |K||w|(Q P+ i)- 
This shows that the (u n ) n and (t>„) n are uniformly bounded in W^^Bp x (0,T)) 
and W^iQp x (0,T)) respectively. 



21 



Step 4. Existence of a solution. 

Now that we know that (u n ) n and (v n ) n are uniformly bounded in compact sets 
with respect to the W 2,1 norm, taking p > N + 1 and using the Morrey inequality, 
we infer that this is also true with respect to the C a norm, for some < a < 1. 
Then, by the Schauder estimates, the time derivative and the space derivatives up 
to order 2 are uniformly Holder continuous in compact sets too. As a consequence, 
((u n ,v n )) n converges (up to subsequences) in C^ c to some (u,v). Passing to the 
limit as n — > oo in (A.l), (A. 2) we eventually find that (u,v) satisfies (3.1)-(3.2). 
Form step 1 we know that u and v are bounded and nonnegative. □ 



B The equation h L (c,f3) = 

In this section we describe in detail how, for £ > small enough, equation 
(6.4) admits two solutions close to r + and r_ respectively. We recall that t± = 
rkiy/ (e/a)£ + 0(£) are the roots of the trinomial g(r) := ar 2 + d^r + e£. Let us 
focus on r+, the other case being analogous. 

Let B be the ball of radius A£, centred at r + ; A large and to be adjusted. For 
r G dB, we have 

\ g (r)\ = a \r-T + \ |r-r_| > aA£(\r + - r_ | - A£) = 2aA^~e~Ja~£ i/2 + 0(£ 2 ). 

On the other hand we have \(p(r, £)| < C£ 3 / 2 + 0(£ 2 ). We can therefore choose A 
large enough and then £ small enough in such a way that \g\ > \tp\ on dB. Since 
g and are holomorphic, by Rouche's theorem the equation (6.4) has the same 
number of solutions in B as g = 0, that is 1. Notice that such a solution has 
positive imaginary part proportional to a/£ and real part of order £. 
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